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ABSTRACT 

Recent investigations have shown that many parameters and assumptions made in the 
application of spectral ageing models to FR-II radio galaxies (e.g. injection index, uniform 
magnetic field, non-negligible cross-lobe age variations) may not be as reliable as previously 
thought. In this paper we use new VLA observations, which allow spectral curvature at GHz 
frequencies to be determined in much greater detail than has previously been possible, to 
investigate two cluster-centre radio galaxies, 3C438 and 3C28. We find that for both sources 
the injection index is much steeper than the values traditionally assumed, consistent with our 
previous findings. We suggest that the Tribble model of spectral ageing provides the most 
convincing description when both goodness-of-fit and physically plausibility are considered, 
but show that even with greatly improved coverage at GHz frequencies, a disparity exists in 
cluster-centre FR-IIs when spectral ages are compared to those determined from a dynamical 
viewpoint. We find for 3C438 that although the observations indicate the lobes are expanding, 
its energetics suggest that the radiating particles and magnetic field at equipartition cannot 
provide the necessary pressure to support the lobes, similar to other cluster-centre source 
such as Cygnus A. We confirm that small scale, cross-lobe age variations are likely to be 
common in FR-II sources and should be properly accounted for when undertaking spectral 
ageing studies. Contrary to the assumption of some previous studies, we also show that 3C28 
is an FR-II (rather than FR-I) source, and suggest that it is most likely a relic system with the 
central engine being turned off between 6 and 9 Myrs ago. 
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1 INTRODUCTION 
1.1 Radio galaxies 

Powerful iFanaroff&Ril^ ( Il974h class I (FR-I) and II (FR-II) 


radio galaxies can have a significant impact on the environ¬ 
ment in which they reside, with structure t hat can extend from 
tens of kiloparsecs jBirkinshaw et al.L |l98ll : [Alexander & Leadi^ 


Il98l Konaretall 

20061; Machalski et all 2009tl to meeanarsecs 

I Mullin et ah, 2006 

; Marscher & Jorsta4,|2011) in size. The abil- 


ity of these radio loud active galaxies to provide the required en¬ 
ergy input to suppress sta r formation in models of galaxy evolution 
through AGN feedback jCroton et all 1 20061 : iBower et"^ l2006ll 


has therefore proved popular in explaining the observed properties 
of stellar populations in massive ellipticals. However, many aspects 
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of the underlying dynamics and energetics of these powerful radio 
sources remain a mystery. 

FR-II radio galaxies generally consist of three large scale 
structures: jets, lobes and hotspots. The jets of FR-IIs are rela¬ 
tivistic and are thought to be the mechanism for transporting ma¬ 
terial from the central active galactic nucleus (AGN) to a termina¬ 
tion shock, usually located at the extremities of the source. This 
shock forms a compact region of synchrotron emission, known 
as the hotspot, and is generally thought to be the dominant loca¬ 
tion of particle acceleration wi thin these sources. As the hotspot 
move s away from the nucleus llBurc 5 1 1971 [I?^: IWinter et all 
ligsd : iMeisenheimer et all [l989l) . the previously accelerated par¬ 
ticles are left behind (possibly with some back flow of mate¬ 
rial) but continue to radiate via the synchrotron process which 
are then observed as the lobes which give FR-IIs their character- 


istic morphology (e. 

e. Scheueil. 1974 

Beaelman & Cioffil 1 19891; 

iKaiser & Alexander!. 1991 Krause et al 

,|2012|). As these lobes are 


in direct contact with the external medium, if we are to determine 
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Table 1. List of target sources and galaxy properties 


Name 

lAU Name 

Redshift 

178 MHz 

5 GHz Core 

Spectral Index 

LAS 

Size 




Flux (Jy) 

Flux (mJy) 

(178 to 750 MHz) 

(arcsec) 

(kpc) 

3C438 

J2153-I-377 

0.290 

48.7 

7.1 

0.88 

22.6 

98.3 

3C28 

J0053-I-261 

0.195 

17.8 

< 0.2 

1.06 

45.6 

148 


‘Name’ and ‘lAU Name’ list the 3C and lAU names of the galaxies as discussed in this paper. ‘Spectral 
Index’ lists the low frequency spectral index between 178 to 750 MHz, ‘LAS’ the largest angular size of the 
source and ‘Size’ its lai'gest physical size. The ‘Redshift’, ‘178 MHz Flux’, ‘5 GHz C ore Flux’, ‘Spectra l 
Index’ , ‘LAS’ and ‘Size’ column values are taken directly from the 3CRR database iLaing et~^1l983h 
(http://3crr.extragalactic.info/cgi/database). 


how these powerful radio sources affect the evolution of galaxies as 
a whole, then we must understand their dynamics, energetics and 
ultimately the total energy they transfer to their environment. 


1.2 Spectral ageing 

In theory, for an electron which is emitting via the synchrotron pro¬ 
cess in a fixed magnetic field, the energy losses scale as 

This leads to a preferential cooling of higher energy electrons and, 
in the absence of any further particle acceleration, produces a spec¬ 
trum which becomes increasingly curved over time. Thus, for elec¬ 
tron energy distribution initially described by a power law 

N{E) = NoE~^ (2) 

we find at later times that 

N{E,9,t) = NoE~\l-E tE)~^~^ (3) 


where Et are the model dependent losses which are a function of 
the pitch angle of the el ectrons to the magneti c field 0, and time 
since acceleration t (see iHarwood et all l2013l for a detailed dis¬ 
cussion). This time dependent process, known as a spectral ageing, 
has become a common method for determining the age of sources 
radiating via the synchrotron process. 

The ability of spectral ageing models to describe the emission 
from the lobes of radio galaxies, particularly their ability to pro¬ 


vide the intrinsic age of a source, has long been a topic of debate 


(e.g. Alexander & Leahvl 1987 ; Eilek & Arendll . [l9^ : lEilek et all 


1 1997 : lBlundell & RawlingsllTOOOh : however, the technique i 
a commonly used tool i n the analysis of both low- and high- 


power radio galax i es (e.g. Jamroz^etal^ 20051: Kharb_et^ , 200^; 
lOrru et all 1 20 id : Ide Gasoerin et all l2012l : iHeesen et al.L 2014h . 
Work to develop more advanced ageing models that bet t er de¬ 


scribe the observed emission continues to this day llTribble , 1^93 


Komissarov & GubanovLll99^:lMachalski et alll2009l:lHardcastle. 


2013h but testing of spectral ageing models against the new gener¬ 
ation of radio interferometers that allow us to obtain much tighter 

constraints remains large ly unexplored. _ 

Our previous work iHarwood et al.L 12013h has given the first 
insights into the spectral structure of FR-II radio galaxies on 
small spatial scales and provided the methods required to allow 
the detailed spectral study of these sources now that truly broad- 
bandwidth observations are available. Whilst this investigation 
found that spectral ageing has at least some basis in reality for pow¬ 
erful radio sources, it has also become apparent that many of the 
previously held assumptions made in the application of models of 


spectral ageing are, at least in some cases, less reliable than previ¬ 
ously thought. 

The first of these potentially incorrect assumptions was that 
the model parameter which describes the observed spectrum of 
the initial power law electron energy distribution (the injection in¬ 
dex) is much steeper than previously assumed. This injection index, 
which is directly related to S by 


S-1 


(4) 


is traditionally assumed to h ave values of around 0 . 5 to 0 .6 based 
on theoretical arguments ( e.g. Blandford & Ostrikeilll978h and ob¬ 
serva tions of hotspots (e.g. lMeisenheimer et al. , 19891 : Carilli et dl 
Il99lh . However, our previous investigation suggests that, for the 
two sources studied, the injection index is >0.8 in both cases. The 
increased low-frequency energy content that this implies has im¬ 
plications for both the spectral ages of powerful radio galaxies and 
their energetics, so it is key that we determine if these findings are 
robust and if they are common to the F R-II populatio n as a whole. 

The second issue, first discussed bv lEilekI jl996h . was a dispar¬ 
ity between the spectral ages and those determined from a dynami¬ 
cal view point. One potential solution to this problem was to place 
tighter constraints on the curvature of the lobe spectrum; however, 
our study showed that even with greatly improved sampling in fre¬ 
quency space, this disparity still remains. In order to resolve this 
issue, it is important to determine whether this age difference is 
found in all galaxies, or if it is confined to FR-IIs in certain envi¬ 
ronments. 

The final problem directly related to the work described within 
this paper is determining which model of spectral ageing best de¬ 
scribes the spectrum of the sources, both in terms of goodness-of-fit 
and in their physical interpretation. We discuss the various models 
of spe ctral ageing further in Section [2.4.2l falso see lHarwood et all 
l2013h but note here that the most commonly a pplied model of spec¬ 
tral ageing proposed bv ijaffe & Perolal lll973l) (the IP model) is fre¬ 
quentl y in conflict with that of KardashevI i ll 9621) and IPacholcz^ 
( Il970l) (the KP model). The JP model is often prefen'ed due to 
its physical plausibility but often provides a worse goodness-of- 
fit than the less physically realistic KP model, and is what we 
also find in our 2013 study. We therefore tested the more recent 
Tribble model of spec tral ageing l lTribblel [19^: iHardcastlel |20 1 3l: 
IHarwood et al.L[^13h which accounts for a more complex descrip¬ 
tion of the magnetic field within the lobes. We found that this model 
was both able to provide a comparable goodness-of-fit to the KP 
model and to retain the physical plausibility of the JP model. If one 
is to determine the impact of powerful radio galaxies on their en¬ 
vironment then resolving these outstanding issues and testing these 
new models in a range of environments in order to be able to re- 
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Figure 1. Combined A, B and C configuration radio map of 3C438 (left) and B and C configuration radio map of 3C28 (right) between 4 and 8 GHz. Imaged 
using multiscale CLEAN and CASA nterms = 2 to a central frequency of 6 GHz (see Section IT^ for details). The off-source RMS of the combined maps is 
8/rJy beam~^ and 4^Jy beam“^ for 3C438 and 3C28 respectively. The restoring beams are 0.30 and 0.99 arcsec for 3C438 and 3C28 respectively and are 
indicated in the bottom left comer of each image. 


liably determine a sources age, hence total power output over its 
lifetime, is a vital step. 


1.3 Outstanding questions 

In order to determine new, reliable parameters and models of spec¬ 
tral ageing, we mus t first investigate whether our previous find¬ 
ings jHarwood et al.L l2013h are intrinsic to a wide range of FR-II 
sources in differing environments, a result of being still too poorly 
constrained in frequency space, or simply outliers in an otherwise 
expected distribution. Within this paper we therefore expand upon 
our initial sample using truly broad bandwidth observations of two 
cluster-centre radio galaxies to address 4 primary questions: 

(i) What are the dynamics and energetics of FR-II sources when 
tightly constrained by broad bandwidth data? 

(ii) Which model of spectral ageing provides the best descrip¬ 
tion of cluster-centre FR-IIs? 

(iii) Is the injection index of cluster-centre FR-IIs flat, as had 
previously been assumed, or steep? 

(iv) Does the well known dynamical versus spectral age dispar¬ 
ity still exist in cluster-core FR-IIs? 

Section [2] gives details of target selection, data reduction and 
the spectral analysis methods used. In Section we present our 
results and Section|4]discusses the implications for our current un¬ 
derstanding of FR-II radio galaxies. Throughout this paper, a con¬ 
cordance model in whi ch Hq — 71 km s~^ Mpc“^, ~ 0.27 
and Da = 0.73 is used dSnergel et Iil [2003) . 


2 DATA REDUCTION AND SPECTRAL ANALYSIS 
2.1 Target selection and observations 


As curvature due to spect ral ageing becomes most easily observ¬ 
able a bove a few GHz jAlexander & Leahvl 1 19871: lAlexandei . 
1 19871 : ICarilli et ^ Il99ll : Perlev. Rbser & MeisenheimeT ' 1997 ; 


iHardcastle & LoonevLl200il) . the resolution obtained with the Karl 
G. Jansky Very Large Array (VLA) at C-band frequencies (4.0- 
8.0 GHz) made it the natu ral choice for the requir e d obs ervations. 
Using the 3CRR sample l lLaing. Riley & Longaij . Il983h of well- 
studied FR-II sources as reference, we selected our targets so as to 
ensure that they were both well sampled by the shortest baselines 
at the highest C-band frequencies, while still maintaining a reason¬ 
able resolution at the lowest frequencies. We therefore restricted the 
largest angular size (LAS) of the sources on the sky to LAS < 50 
arcsec at a redshift of z < 0.5. We also imposed a spectral index 
cut-off for the integrated spectrum of q > 0.8, where 0.8 is the av¬ 
erage of the 3CRR sample. This ensures that the targets are likely 
to show significant spectral ageing rather than being very young, 
nearby sources, selected due to the small apparent size of the LAS 
cut. 

From the eight sources that satisfied these criteria in 3CRR, 
we excluded the peculiar sources 3C171 and 3C305, where the 
radio structure plausibly re sults from strong interac tion with cold 
gas within the host galaxy dHardcastle et al.L ItOIOI) . and 4C 14.27 
which is likely to be a relic system. 17.5 hours of observing time 
was awarded which was suitable for observation of 2 of the re¬ 
maining 5 sources in the required configurations. We therefore split 
the time equally between 3C438 and 3C28, both of which reside 
in cluster-centre environments. Details of the general properties of 
these sources are given in Table[T] 

In order to ensure that all known compact and diffuse emission 
was sampled, observations were made using the A, B and C array 
configurations. As 3C28 lacks any compact structure on scales that 
are sampled by the A configuration array, only B and C were used 
in the final images. At the time of observation the maximum band¬ 
width capability of the VLA was 2 GHz as the 3-bit samplers were 
not available; we therefore used two 2 GHz observations to give 
full frequency coverage between 4 and 8 GHz. 


2.2 Data reduction 

The data were downloaded from the NRAO archive and reduced 
using CASA in the standard manner as described in the CASA cook- 
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boo]|3 and the online tutorial for narrow-field, C-band observations 
0. As multiple images were required for our analysis rather than 
just a single broadband data cube, each of the 32 125-MHz spectral 
windows was individually calibrated. 

The structures of 3C438 and 3C28 differ significantly, hence 
different imaging strategies were required to obtain the optimal 
radio maps. For 3C438, the A-configuration data were first self- 
calibrated in phase only and CLEANed to convergence. The A con¬ 
figuration images were then used as a model to cross-calibrate the 
B-configuration observations. The A and B configuration uv data 
were then combined and imaged, and these images in turn were 
used to cross-calihrate the C configuration data. As 3C28 does not 
contain compact structure, it was not possible to apply the same the 
cross-calibration method used for 3C438. Instead, the B and C con¬ 
figuration uv data were combined and a standard self-calibration 
performed in phase only. The radio maps of both sources were pro¬ 
duced using the multi-scale cl ean algorithm MS CLEAN jCornwelH 
|2008 LI^u & Cornwelll20llh at scales of 0, 5, 15 and 45 times the 
cell size, where the cell size was set to one fifth of the beam size. 
In order for a detailed spectral analysis to be performed the param¬ 
eters of each image must be matched, so the beam size of the final 
images was set to that of the lowest frequency data. We excluded 
the two lowest spectral windows for both sources (for reasons dis¬ 
cussed in Section and so image at the resolution of the lowest 
frequency used in the analysis, rather than the lowest frequency ob¬ 
served. A summary of parameters used for the imaging described 
in this section is shown in Table 

As the observations described so far were taken for each array 
configuration in a single pointing we can be confident that they are 
inherently well aligned. However in order to correct for small, sub 
pixel variations in position, Gaussians were fitted to the radio core 
for 3C438 and to a point source close to 3C28 (where the radio 
core is not clearly visible in all images) using CASAVIEWER. A 
reference pixel was then decided upon and each map aligned using 
the AIP^ OGEOM task. 

In order to extend our frequency coverag e the combined VLA / 
MERLIN L-band data of 3C438 presented bv lTreichel et all l l200lh 
and data for 3C28 from the atlas of DRAGNs (Leahy, Bridle & 
StrorrQ), which are of comparable resolution to the C-band obser¬ 
vations, were also obtained. These data were already calibrated to 
a high standard, hence only regridding to J2000 coordinates, re¬ 
imaging to match the parameters of Table and alignment as de¬ 
scribed above was required. 

To check the impact of any residual alignment errors, Gaus¬ 
sians were again fitted to the aligned maps. We find the standard 
deviation in the alignment between maps to be only 0.007 pixels, 
and so misalignments are unlikely to make any significant contribu¬ 
tion to the uncertainty in the fitting of models to the oldest regions 
of plasma. 

Additional images were also created using the full bandwidth 
centred at a frequency of 6 GHz. As, for the combined data set, the 
flux changes significantly over the bandwidth, the data must be suit¬ 
ably scaled in order to produce a realistic image of the source. We 
therefore used the CASA multi-frequency synthes is (MFS) CLEAN 
parameter nterms = 2 l lRau & Comwelj 1201 ih which scales the 


^ http://casa.nrao.edu/docs/UserMan/UserMan.html 

^ http://casaguides.nrao.edu/index.php?title=EVLA_6cmWideband_Tutorial 
.SN2010FZ 

® http://www.aips.nrao.edu/index.shtml 
http://www.jb.man.ac.uk/atlas/ 


Table 2. Summary of imaging parameters 


Source 

Parameter 

CASA Name 

Value 

Units 

3C438 

Polarization 

STOKES 

I 



Image Size 

IMSIZE 

4096 4096 

Pixels 


Weighting 

ROBUST 

-0.5 



Cell Size 

CELL 

0.084 0.084 

Arcsec 


Beam size 

RESTORINGBEAM 

0.42 0.42 

Arc sec 


Multiscale 

MULTISCALE 

[0, 5, 15, 45] 

Pixels 

3C28 

Polarization 

STOKES 

I 



Image Size 

IMSIZE 

6144 6144 

Pixels 


Weighting 

ROBUST 

1 



Cell Size 

CELL 

0.42 0.42 

Arcsec 


Beam size 

RESTORINGBEAM 

2.1 2.1 

Arcsec 


Multiscale 

MULTISCALE 

[0, 5, 15, 45] 

Pixels 


‘Parameter’ refers to the imaging parameter used in making of radio maps 
within this chapter. ‘CASA Name’ refers to the CASA parameter taking the 
value stated in the ‘Values’ column. 


flux by a spectral index value fitted over the observed bandwidth. 
The resulting images are shown in Figure[T] 

2.3 Data quality 

From the 32 VLA spectral windows between 4 and 8 GHz, we 
found that 28 had good quality data in all array configurations for 
both 3C438 and 3C28. Of the remaining four, the spectral win¬ 
dow at 6.10 GHz (SPW 17) was of extremely poor quality in all 
configurations and was therefore excluded from our analysis. At 
4.02 GHz (SPW 0), heavy flagging of RFI was required in all array 
configurations leading to poor data quality and so this image was 
excluded from further analysis. The observations at 4.14 and 6.36 
GHz (SPW 1 and 19 respectively) are of reasonable quality, but sig¬ 
nificant RFI was present in the B and C configuration data. These 
two frequencies were therefore excluded from the spectral analysis, 
but included in the making of the combined image shown in Fig- 
ure[T] As these problems are common to both sources, we suggest 
that they are most likely due to an instrumental issue at the time of 
observation. It is also worth noting that the low flux density ‘hole’ 
observed in the northern lobe of 3C43 8 is also seen in t he atl as of 
DRAGNs and in observations made bv iHardcastle et all ( ll997h and 
so is likely to be physical in origin, rather than an imaging artefact. 

Within the field of view of our observations of 3C28 an¬ 
other large, bright, extended radio source was also serendipitously 
observed. Unfortunately, due to being located at the edge of the 
JVLA’s primary beam and with no well matched L-band data avail¬ 
able, we are unable to perform a full spectral analysis. However, 
we are still able to provide some insight into the source’s proper¬ 
ties and therefore provide a brief analysis of this peculiar source in 
Appendix lAl 


2.4 Spectral analysis 

For the analysis of our sources we have used the Broadband Ra¬ 
dio Astronomy Toolfl (BRATS) software package which provides a 
wide range of spectral analysis and spectral age model fitting tools. 
The basic usage and functionality of BRATS has been discussed in 

® http://www.askanastronomer.co.uk/brats 
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detail previously l lHarwood et all l2013h and in the BRATS cook- 
bool 0 and so we do not provide an in depth discussion here. We 
instead give a brief summary of the fitting performed and high¬ 
light new features which have been added to the software and used 
within this paper over previous works. 

2.4.1 Region selection 

A region fully encompassing the source (but excluding the core) 
was defined using along with a background region for de¬ 

termining the off-source thermal noise. To test the impact of the 
strong emission in 3C438 coincident with the assumed location of 
the jet on the model parameters and fits, a second region was also 
created with this area excluded. Initial source detection was then 
performed, for which a 5cr cut-off was us ed based on thi s RMS 
value. We adopt the same practice used bv iHarwood et alj l l2013tl 
and apply an on-source multiplier of 3 times this value for the deter¬ 
mination of statistical values and region selection to account for the 
increased RMS noise due to the increased uncertainty in the mod¬ 
elling of the extended emission. As each VLA image is taken within 
a single pointing, one can be confident that the relative flux calibra¬ 
tion error between images is small. We therefore take the flux cal¬ 
ibration error to be 1 per cent for the VLA images between 4 and 
8 GHz. As the L-band observations were made with MERLIN and 
the old VLA system, the flux calibration error relative to the VLA is 
not immediately clear. We therefore take a conservative value of the 
standard absolute flux calibration en'or value for VLA L-band ob¬ 
servations of 5 per cent. To reduce the impact of the superpositio n 
of spectra on our results jHarwood et al.Ll2013l : IStroe et al.Ll2014h . 
and as these data are of good quality, in the analysis that follows 
we considered each source on a pixel by pixel basis. A summary of 
the values chosen for the region selection is shown in Table|4] 


Table 3. Summary of image noise by frequency 


Frequency 

(GHz) 

3C438 

Off-source On-source 

(^Jy/beam) (/iJy/beam) 

3C28 

Off-source On-source 

(/iJy/beam) (/iJy/beam) 

1.42 

70.6 

212 

- 

- 

1.46 

- 

- 

75.3 

226 

4.33 

34.1 

102 

16.6 

49.8 

4.45 

33.5 

101 

15.8 

47.4 

4.58 

35.6 

107 

16.3 

48.9 

4.71 

34.1 

102 

16.2 

48.6 

4.84 

33.8 

101 

14.4 

43.2 

4.97 

32.1 

96.4 

14.4 

43.2 

5.07 

38.8 

116 

15.5 

46.5 

5.20 

33.5 

100 

14.2 

42.6 

5.33 

32.2 

101 

14.7 

44.1 

5.45 

30.9 

96.7 

15.1 

45.3 

5.58 

28.8 

92.7 

15.0 

45.0 

5.71 

28.7 

86.4 

14.6 

43.8 

5.84 

28.7 

86.0 

15.4 

46.2 

5.97 

31.6 

94.8 

16.2 

48.6 

6.03 

36.0 

108 

11.9 

35.7 

6.29 

25.8 

77.3 

9.77 

29.3 

6.54 

25.7 

77.0 

9.63 

28.9 

6.67 

27.6 

82.9 

10.6 

31.8 

6.80 

25.5 

77.6 

9.39 

28.2 

6.93 

26.3 

79.0 

9.40 

28.2 

7.03 

29.6 

88.8 

10.2 

30.6 

7.16 

23.9 

71.7 

9.12 

27.4 

7.29 

23.8 

71.3 

9.31 

27.9 

7.41 

23.9 

71.8 

9.52 

28.6 

7.54 

25.7 

77.1 

9.85 

29.6 

7.67 

24.8 

74.5 

9.54 

28.6 

7.80 

23.8 

71.4 

9.76 

29.3 

7.93 

24.7 

74.2 

9.45 

28.4 


2.4.2 Model fitting and parameter determination 


We have tested the broad bandwidth data described above against 
three single-injection models of spectral ageing. The two most 
common models used to describe emission from the lobes of 

and 


FR-II galaxies ar e those proposed by [Kardashe^_([l9ffi ) 
IPacholcz^ ( Il97(lh (the KP model) and by Jaffe & Perolal ' 




(the JP model). Thes e models, which have been wide l y dis- 
cussed elsewhe re (e. g. iLeah^ Il99ll : ISlundell & RawlingsLl 200 l] : 
iHarwood et alil^l3h . differ only in their treatment of the pitch an¬ 
gle of electrons with respect to the magnetic field. The KP model 
assumes the angle of any given electron is fixed over its radiative 
lifetime, compared to the JP model which uses a time averaged 
value. This time averaging leads to an exponential cut-off at high 
frequencies, compared to the KP model which is comparatively flat 
due to there being a supply of high energy electrons at small pitch 

angles which are capable of radiating at higher fre quencies. _ 

The third model tested is that proposed bv lTribbl^ jl993tl 
which attempts to account for a spatially non-uniform magnetic 
field. Both the KP and JP models assume a fixed magnetic field 
strength but this is unlikely to be physically realistic. The Tribble 
models therefore attempt to account for the magnetic field structure 
by assuming it to be Gaussian random and all owing elect r ons to 
diffuse across regions of varying field strength. lHardcastI3 l l2013tl 
has recently shown that in the weak field, high diffusion case where 


® http://www.askanastronomer.co.uk/brats/bratscookbook.pdf 
^ http://ds9.si.edu 


‘Frequency’ refers to the frequency of the map, ‘Off-source RMS’ refers to 
RMS noise measured over a large region well away from the source and 
‘On-source RMS’ the noise used for region selection and statistics as per 
Section im 


electrons are free-streaming, the spectrum can be modelled by inte¬ 
grating the standard JP losses over a Maxwell-Boltzmann distribu¬ 
tion in field energy density. We have previo usly detailed the appli- 
cation of this model to similar FR-II sources jHarwood et allliof 3h 
suggesting that it may provide a solution to an outstanding problem 
where in many cases the KP model provides a better goodness-of- 
fit than the more physically realistic JP model. 

We have made the standard assumption that the magnetic 
field (or in the case of the Tribble model, mean magnetic field) 
is in equipartition, which we calculate u sing the SYNCH code of 
iHardcastle. Birkinshaw & Worralll (Il998h . We took the minimum 
and maximum electron Lorentz factor to be 7 = 10 and 7 = 
1 X 10® respectively. To determine the best fitting injection index 
we have used BRATS’ findinjecf command between ainj — 0.5 
and 1.0 for 3C438 and between 0.5 and 1.5 for 3C28 at intervals of 
0.1. A second run was then performed around the minimum with 
a smaller step size of 0.01. The resulting values were then plotted 
and the minimum injection index found. Final model fitting of the 
sources was then performed in BRATS using the derived values. The 
statistical values for each model were recorded and spectral ageing, 
and error maps were then exported ready for further analysis. 
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Table 4. Summary of adaptive region parameters 


Parameter 

Value 

Units 

Description 

Signal to noise 

I 


SNR (pixel to pixel) 

Search area 

1 

Pixels^ 

Max. search area 

On-source noise 

3 


On-source noise multiplier 

Hot pixel limit 

20 

Per cent 

Max. pixel vaiiation 

Map variations 

-1 


Maximum map variation (off) 


‘Value’ refers to the values applied within BRATS for the corresponding ‘Pa¬ 
rameter’. The ‘Description’ column provides further details of the value 
meaning. Note that a signal to noise of 1 is still subject to the Scr cutoff 
and the on-source noise multiplier. 

2.4.3 Uncertainty calculation 

One key feature which has been applied to the data over previous 
investigations is the inclusion of improved uncertainty calculations. 
During the fitting process BRATS performs a search over a large 
range of potential spectral a ges allowing the curve as a func¬ 
tion of age to be determined. lAvnil ( Il976h shows that the Icr error 
of such a minimization is given by a deviation of Ax^ ~ 1 from 
the minimum x^ value. These errors were therefore calculated dur¬ 
ing the fitting process, mapped as a function of position, and then 
exported for further analysis. A similar error calculation was also 
performed to provide Icr errors for the injection index values. How¬ 
ever, as here one uses the sum of x^ over all regions for each injec¬ 
tion index to determine the minimum value, they are over weighted 
by a factor of the beam area. Assuming the injection index is ap¬ 
proximately constant across any given source, the injection index 
errors are therefore determined by finding where Ax^ = 1 for the 
corrected x^ values such that xlorr = I-beam, where Ab^am is 
the beam area. These methods provide a much more robust error 
estimate than in previous studies and are applied throughout this 
paper. 

2.4.4 Spectral index fitting 

An additional improvement to BRATS is the enhanced handling 
of spectral index fits. As the uncertainty on the flux measure¬ 
ments of radio observations varies between images, observations 
and telescopes, the assumption of constant errors used by a standard 
least-squares fit is often not valid. BRATS therefore now includes 
a weighted least-squares option for spectral index determination. 
This function fits a linear regression in log space of the standard 
form y(c, x) = cq -\- c\x but, instead of assumiM uniform errors, 
uses the GSL ‘wlinear’ and ‘linear_est’ function^ to minimize the 
weighted sum of squared residuals such that when fitting to n ob¬ 
servations 

n 

X^ = ^ Wi{yi - (cq + ciXi)f (5) 

i=l 

where the weights are given by w = 1/crf and Oi is the error on a 
given flux measurement. Note that as fitting is performed in log-log 
space, the fractional errors on the measurements must be small so 
that one can assume they remain approximately symmetric when 
transformed to natural log space. In general, this assumption is 

® https://www.gnu.org/software/gsl/manual/htmLnode/Linear- 
regression.html#Linear-regression 


Table S. Best fitting injection indices and magnetic field strengths 


Source 

Model 

Injection 

Index 

En'or 

-1- 

Magnetic Field 
Strength (iiT) 

3C438 

JP 

0.84 

0.01 

0.01 

4.02 


KP 

0.78 

0.01 

0.01 

3.65 


Tribble 

0.82 

0.01 

0.01 

3.89 

3C438 

JP 

0.80 

0.01 

0.01 

3.77 

(No Jet) 

KP 

0.72 

0.01 

0.02 

3.34 


Tribble 

0.77 

0.01 

0.01 

3.60 

3C28 

JP 

1.21 

0.01 

0.01 

1.35 


KP 

1.12 

0.02 

0.02 

1.12 


Tribble 

1.17 

0.02 

0.01 

1.22 


Best fitting injection indices for 3C438 and 3C28 and magnetic field 
strengths. ‘No Jet’ indicates the best fitting injection index values where 
emission coincident with the assumed location of the jet is ex cluded (Sec- 
tion l2.4.l1 and l3.U . EiTors are determined using the methods of lAvnil jl97fSh 
detailed in Section l2.4.3l 

valid for fractional errors <10 per cent and so is well within the 
uncertainties used within this paper. 

3 RESULTS 

3.1 Model Parameters 

The plot of x^ values for varying injection indices are shown 
in Figure 1^ where one can see that for both sources the minima 
are model-dependent and occur over a range of values (Table |5](- 
This is in contrast to our earlier results on other FR-II sources 
jHarwood et al.L 12013h where we used the same method of deter¬ 
mining the injection index but found that all models minimized to 
the same value. The most likely cause of this difference is the use 
of truly broad bandwidth observations, providing 29 data points 
ov er the frequency spac e compared to the maximum of 9 used 
bv iHarwood et alj ( l2013h . As the curvature at GHz frequencies is 
more tightly constrained, one is able to differentiate between much 
smaller variations in the models which in turn leads to a variation in 
the derived injection index values; the observed model dependence 
is thus not entirely unexpected. 

Regardless of the model used, one thing which is immediately 
clear is that the injection index of both sources remains steeper than 
the traditionally assumed values of 0.5 to 0.6. This is particularly 
prominent in 3C28 where the injection index is found to be between 
1 .12 ±o;o2 (KP) and 1.21±o;oi (JP), up to double the previously as¬ 
sumed values. Although not as dramatic as those seen in 3C28, the 
injection index of 3C438 is also found to be relatively steep with 
values ranging between 0.78±o;o2 (KP) and 0.84±o;oi (JP) when 
the whole source is co nsidered. However, un like 3C28 or those 
sources investigated by iHarwood et al.l l l2013h . 3C438 has bright 
emission coincident with the assumed location of the jet (Figure[T]( 
which complicates the determination of the injection index. 

The jets of FR-IIs are still poorly understood and it is there¬ 
fore not clear whether the plasma responsible for this emission is 
due to reaccleration of lobe material at a jet-lobe boundary, or if it 
is accelerated internally within the jet itself. In either case, there is 
no a priori reason to believe that the observed jet emission should 
be well described by a spectral ageing model. We therefore ran 
a second minimization with the jet excluded from consideration. 
From Figurel^and Table|5] one can clearly see that with the jet ex¬ 
cluded a minimum occurs at a decreased injection index of between 
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Figure 2. values for 3C438 including the jet emission (top left), excluding the jet emission (top right), and 3C28 (bottom) for varying injection index values 
using the JP, KP and Tribble models of spectral ageing. Data points are taken at intervals of 0.1 between 0.5 and 1.0 for 3C438 and 0.5 and 1.5 for 3C28, with 
intervals of 0.01 around the minimum values. The data are fitted with natural cubic splines. As all points lie on the fitted spline they are excluded for clarity. 
The best fitting injection index values are shown in Table|^ 


0.72±o!o2 (KP) and 0.80±o;oj (JP). At first glance, this reduction 
in injection index is counter-intuitive as these bright regions are 
thought to be due to freshly accelerated plasma and so should con¬ 
tain the flattest spectrum emission. However, we see from Figure|3 
that a large fraction of the emission has a spectral index of q > 1 . 
As these jet regions are also described by the ageing models (even 
if incorrectly) as having a low spectral age (see Section 13.21 and 
Figure |4} they are not subject to additional curvature at low fre¬ 
quencies and maintain their steep, power law spectra driving the 
injection index to higher values. Given the non-negligible effect of 
the jet on the injection index, we perform all subsequent analysis 
with the jet regions both included and excluded. 


For both sources, we use the best fitting injection index val¬ 
ues of Table |5] for the final model fitting and determination of the 
equipartition magnetic field strength as detailed in Section [2.4.21 
with an initial electron energy power law index of 5 = 2onni + 1 . 
We find a magnetic field strength in the range of 3.65 nT (KP) to 
4.02 nT (JP) for 3C438, 3.34 nT (KP) to 3.77 nT (JP) for 3C438 
where the jet has been excluded and 1.12 nT (KP) to 1.35 nT (JP) 
for 3C28. These magnetic field strengths, summarised in Tabled 
are used for all subsequent analysis. 


3.2 Spectral age and model comparison 

3.2.1 3C438 

Figures |4] and |5] show the spectral age, and error values as a 
function of position for the JP, KP and Tribble models for 3C438. 
We see that, even when the obvious variations due to the jet are ig¬ 
nored, significant cross-lo be age variations exist in a manner con¬ 
sistent with the findings of lHarwood et al.l ( 12013l) . Low age regions 
in the northern and southern lobes are coincident with the hotspots 
as one would expect, with older regions of plasma residing closer 
to the core and the main spectral features remaining independent 
of model type. With the jet regions excluded (Figure|5} we see that 
the spatial distribution of ages mirrors that of the full source as ex¬ 
pected. From Table |7] one sees that for all models the spectral ages 
have increased in the case where the jet is not considered due to the 
steeper injection indices and lower magnetic field strengths. 

A particularly interesting spectral feature which is observed 
in Figures |4] and [5] is the low age region at the edge of the south¬ 
ern lobe. It is not immediately obvious why particle acceleration 
should be occurring over such an extended region of the source, so 
further testing is requ ired. We therefore use d BRATS and the L- and 
X-band presented bv lTreichel et al.l l l200lh . reimaged at 0.42 arc- 
sec resolution, to produce a two point spectral index map. Figure| 6 ] 
shows the resultant image zoomed to the southern lobe. Although 
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Table 6. Model Fitting Results 


Source 

Model 

Mean 

Mean , 

'^red 

< 68 

Confidence Bins 

68-90 90-95 95 -99 

> 99 

Rejected 

Median 

Confidence 

3C438 

JP 

47.17 

1.75 

2229 

1214 

498 

1203 

4108 

Yes 

> 99 


KP 

45.87 

1.70 

2416 

1228 

546 

1175 

3887 

Yes 

> 99 


Tribble 

46.85 

1.74 

2256 

1237 

513 

1196 

4050 

Yes 

> 99 

3C438 

JP 

45.05 

1.67 

1779 

1019 

449 

956 

2930 

Yes 

> 99 

(No Jet) 

KP 

43.16 

1.60 

2013 

1053 

472 

932 

2663 

Yes 

> 99 


Tribble 

44.42 

1.65 

1828 

1053 

448 

951 

2853 

Yes 

> 99 

3C28 

JP 

34.45 

1.42 

1688 

307 

126 

164 

644 

No 

< 68 


KP 

35.04 

1.17 

1711 

280 

90 

164 

684 

No 

< 68 


Tribble 

34.41 

1.27 

1705 

300 

108 

162 

654 

No 

< 68 


‘Model’ refers to the spectral ageing model fitted to the target listed in the ‘Source’ column. Mean lists the average over the entire source with an 
equivalent reduced value shown in the ‘Mean X^ed column. ‘Confidence Bins’ lists the number of regions for which their y^ values falls with the stated 
confidence range. ‘Rejected’ lists whether the goodness-of-fit to the source as a whole can be rejected and ‘Median Confidence’ the confidence level at which 
the model can or cannot be rejected. 


St^ectral Index I 1.5 2 



’=d- 

CsJ 


o 

L_I_I_I_I_I_I_ \ _^_I 

£S'’55'"52?B 52?5 52?4 52f3 52=2 52?1 52?Q 51?9 

r; 

Right Ascension 

Figure 3. Spectral index map of emission coincident with the assumed loca¬ 
tion of the jet in 3C438. The spectral index fitting is performed between 1.4 
and 7.9 GHz using the weighted least-squares method described in Section 

a small region of low spectral index is observed, there is a distinct 
lack of the extended fiat spectrum emission which one would ex¬ 
pect if significant particle acceleration was occurring. The maps 
also provide some insight into this spectral feature. The small area 
of flat spectrum / low age seen in both the spectral index and age¬ 
ing maps is a very poor fit, with the remainder of this region having 
higher values compared to the source as a whole. From Figure 
|T]one can see a sharp gradient in flux on which this edge lies and 
it is likely that this low age region is due to the sharp edge of the 
lobe. We therefore consider the spectral ages derived not to repre¬ 
sent the true age of the plasma within this elongated region for the 
remainder of this paper. 

In addition to the lobes, 3C438 has significant emission co¬ 


incident with the assumed location of the jet. As was discussed in 
Section [TT] there is no reason to believe that it should be well de¬ 
scribed by a spectral ageing model; however, it does provide a good 
proxy for changes in the steepness of this emission’s spectrum. In 
the southern lobe, changes in the spectrum are observed to be rel¬ 
atively smooth with an increase in brightness and flattening of the 
spectrum occurring as the jet passes into the observed lobe emis¬ 
sion. At this same location, a change in direction of the jet emission 
is also seen to occur. In the northern lobe, the observed jet emission 
remains straight along its length, but is seen to be ‘knotty’ with 
sharp changes in the spectrum. The errors associated with these re¬ 
gions (Figure |4ll are higher than those for the surrounding plasma; 
however, from the maps one can see that the bright regions of 
jet emission are well fitted by spectral ageing models, but are a 
poor fit in lower flux regions. This is perhaps not surprising as, if 
the accelerated particles take the form of a power law, it will be 
well described by a spectral ageing model of zero age, regardless 
of any subsequent processes. We discuss this jet emission further 
in Section lA^ 


Although the distribution of ages is consistent with what one 
would expect from spectral ageing, the goodness-of-fit of the mod¬ 
els tested is far less clear. From Figure|4]one can see that the overall 
y^ values are fairly consistent over the lobe emission ranging be¬ 
tween Xreducsd ~ 1 ftid 5 (excluding the southern edge described 
above). The northern lobe provides a noticeably better fit compared 
to the southern lobe, which may give a hint towards the cause of 
these high values; we discuss this point further in Section 
The northern hotspot is also poorly fitted by all 3 models of spec¬ 
tral agein g, a feature which was also common in the sources inves¬ 
tigated bv iHarwood et al.l ( l2013h . The two southern hotspots do not 
show these high values of but they are much weaker than their 
northern counterpart. From Table [ 6 ] we find that for the goodness- 
of-fit of the 3 models, KP provides the best fitting model, the IP 
model providing a noticeably poorer model of the source and the 
Tribble model providing an intermediate goodness-of-fit. However, 
all 3 models are rejected at the 99 per cent significance, which given 
the overall consistency of the fit to the lobes is a somewhat surpris¬ 
ing result. We discuss the possible cause of these systematically 
high y^ values further in Section lA^ 
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Figure 4. Spectral ageing maps (left), corresponding maps (middle) and error maps (right) of 3C438 with 7.9 GHz flux contours. Three model fits are 
shown; JP model (top), KP model (middle) and Tribble model (bottom) using the best fitting injection index given in Table[5] 
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Table 7. Lobe Advance Speeds 


Source 

Model 

Lobe 

Max Age 
(Myrs) 

+ 

- 

Distance 

(kpc) 

Speed 
(10~^ v/c) 

+ 

- 

3C438 

JP 

North 

3.00 

0.05 

0.08 

30.45 

3.313 

0.052 

0.088 


KP 

North 

4.35 

0.40 

0.29 

30.45 

2.285 

0.210 

0.152 


Tribble 

North 

3.61 

0.10 

0.07 

30.45 

2.753 

0.076 

0.053 


JP 

South 

3.00 

0.05 

0.08 

25.45 

2.769 

0.046 

0.074 


KP 

South 

4.35 

0.40 

0.29 

25.45 

1.909 

0.176 

0.127 


Tribble 

South 

3.61 

0.10 

0.07 

25.45 

2.301 

0.064 

0.045 

3C438 

JP 

North 

3.36 

0.07 

0.07 

30.45 

2.958 

0.062 

0.062 

(No Jet) 

KP 

North 

5.50 

0.95 

0.41 

30.45 

1.807 

0.312 

0.135 


Tribble 

North 

4.19 

0.10 

0.10 

30.45 

2.372 

0.057 

0.057 


JP 

South 

3.36 

0.07 

0.07 

25.45 

2.472 

0.052 

0.052 


KP 

South 

5.50 

0.95 

0.41 

25.45 

1.510 

0.261 

0.113 


Tribble 

South 

4.19 

0.10 

0.10 

25.45 

1.982 

0.047 

0.047 

3C28 

JP 

North 

12.42 

0.44 

0.44 

59.63 

1.567 

0.056 

0.056 


KP 

North 

18.38 

1.10 

0.84 

59.63 

1.059 

0.063 

0.048 


Tribble 

North 

16.82 

0.65 

0.64 

59.63 

1.157 

0.045 

0.044 


JP 

South 

12.42 

0.44 

0.44 

62.86 

1.652 

0.059 

0.059 


KP 

South 

18.38 

1.10 

0.84 

62.86 

1.116 

0.067 

0.051 


Tribble 

South 

16.82 

0.65 

0.64 

62.86 

1.220 

0.047 

0.046 


‘Model’ is the spectral ageing model fitted to the target listed in the ‘Source’ column. ‘Max Age’ is the maximum age of the corresponding ‘Lobe’ in Myrs. 
‘Distance’ gives the separation between the oldest aged population in the lobe and the hotspot in kpc (note this is not necessarily the length of the lobe itself). 
‘Speed’ lists the derived advance speed as a fraction of the speed of light as detailed in Section [T3l Note that the lobe advance speeds of 3C28 contain additional 
significant uncertainties as detailed in Sections [3.3l andl^ 


3.2.2 3C28 

Figure |7] shows the spectral age, ^ and error values as a func¬ 
tion of position for the JP, KP and Tribble models for 3C28. While 
the general distribution of ages is what one would expect for an 
FR-II source (i.e. a gradient in age from the tip to the core). We 
again see two zero age regions located at the edge of the source 
where no obvious source of particle acceleration is present. These 
non-physical, low age regions have now been observed in all FR-II 
galaxies since the method’s development and are therefore likely 
an unavoidable consequence of fitting on such small spatial scales 
in such high dynamic range sources. However, as we are primarily 
concerned in this paper with the general spatial distribution of ages, 
overall goodness-of-fit, and the oldest regions of plasma, these non¬ 
physical regions have a minimal impact on the analysis and conclu¬ 
sion drawn and so are excluded from subsequent analysis. 

Comparing Figure |7] to the spatial distribution of ages in 
3C438 it is evident that the gradient of ages is far smoother, lacking 
any compact structure or jet interaction. There are also no low age 
regions normally indicative of particle acceleration, with youngest 
age region of the source on the order of between 6 (JP) and 9 (KP) 
Myrs with a maximum age of approximately 12 (JP) and 18 (KP) 
Myrs. Comparing the spectral ages to the flux map of Figure [T] 
this is perhaps not too surprising for the northern lobe where no 
compact structure is observed; however, the southern lobe contains 
what has so far been assumed to be an active (albeit weak) hotspot. 
This lack of obvious low age, particle acceleration regions is a sig- 
ni ficant difference comp ared to 3C438 and those sources studied 
by iHarwood et alj ( 1201 3h and may have implications for our cur¬ 
rent understanding of 3C28’s current position in its radio life-cycle. 
We discuss this further in Section l43] 

From the maps of Figure [7] we see that 3C28 is well fit¬ 
ted across the majority of the source with reduced between 1 
and 2. One particularly interesting feature is the high regions 
in the northern lobe peaking at Xreduced ~ 8> similar to the high 


X^ values we f ound for the hotspots of 3C300 and 3C436 in our 
previous work dHarwood et al.L l2013l) . From Table we find that 
the goodness-of-fit of the 3 models follows the same sequence as 
3C438 with the KP model providing the best fit, the JP model pro¬ 
viding a noticeably poorer fit and the Tribble model providing an 
intermediate goodness-of-fit. However, in contrast to the results for 
3C438, all of the models provide a good description of the observed 
spectrum as a whole with none of the 3 models being rejected at 
even the 68 per cent confidence level. 


3.3 Lobe Speeds 


To determine the characteristic ad vance speed of the lo b es we use 
the standard method employed by Alexander & Leah^ ( Il987h . as 
in our previous study dHairwood et aU l2013l) . Assuming an electron 
population accelerated at the hotspot, the advance speed of the lobe 
can be given by viobe = dha/toU, where viobe is the characteristic 
lobe speed, dhs is the distance to the current location of the hotspot 
and told is the age of the oldest region of plasma. Note that for 
the sources presented within this paper we use the term hotspot to 
refer to the regions of bright emission at the end of the radio lobes, 
although we acknowled ge that this does no t fit the strict definition 
laid out by, for example. iLeahv et al.l(ll997h . 


The advance speeds for both sources are shown in Table |7] 
For 3C438 and using the best fitting injection index we find an ad¬ 
vance speed of between 2.29 (KP) and 3.31 (JP) xl0“^ c for the 
northern lobe and 1.90 (KP) and 2.77 (JP) x 10“^ c for the south¬ 
ern lobe. These values decrease slightly where the jet is excluded to 
between 1.81 (KP) and 2.96 (JP) x 10“^ c for the northern lobe and 
1.51 (KP) and 2.47 (JP) x 10“^ c for the southern lobe. These ages 
fall well within the expected range compared to previous studies of 
similar sources (e.g. [Alexander & LeahvLll987h . 

For 3C28 we derive lobe advance speeds slower than those 
found in 3C438 with values between 1.06 (KP) and 1.57 (JP) 
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Figure 5. Spectral ageing map (left), con'esponding map (middle) and error map (right) of 3C438 with emission coincident with the assumed location of 
the jet excluded and overlaid with 7.9 GHz flux contours. As all models follow the same spatial distribution of ages shown in Figure|4] only the JP model ht is 
shown here for conciseness. The best htting injection index given in Table|^is used with the maximum spectral ages of all models given in Table|7] 


X10“^ c for the northern lobe and 1.12 (KP) and 1.65 (JP) x 10“^ 
c for the southern lobe. Whilst these values appear reasonable com¬ 
pared to previous studies of FR-II sources, the method used may in 
this case be significantly underestimating the lobe advance speed. 
The lack of compact emission and, perhaps most significantly, low 
age emission in 3C28 mean that both dhs and told are uncertain. In 
the northern lobe, the distinct lack of a hotspot means that deter¬ 
mining the exact location of the initial particle acceleration is dif¬ 
ficult; however, given some compact structure typical of a hotspot 
is observed in the southern lobe (Figure [TJ the assumption that the 
sites of particle acceleration in 3C28 are located near the tip of the 
lobes (determined here by the peak flux), as is observed in the ma¬ 
jority of FR-II sources, is not unreasonable. Therefore whilst the 
value of dhs is uncertain, it is unlikely to change the advance speed 
by a large amount. 

The validity of told may have a much larger impact on the de¬ 
rived speeds of 3C28. All of the characteristic advance speeds con- 
sidered so far (3C438. [Alexander & Leah^ 1 19871 iHarwood et all 
l2013h have been for sources which are cun'ently active and advanc¬ 
ing through the external medium. However, the lack of observed 
low age emission (Figure |7]( means that this may not be the case 
for 3C28. Under the assumption that the source is no longer ac¬ 
tive, the time variable for determining the advance speed becomes 
told = tmax — tmin whete tmax and tmin ate the maximum and 
minimum observed spectral ages respectively, increasing the char¬ 
acteristic advance speed of 3C28 during the active phase by a factor 
of around 2. We discuss the possibility of 3C28 being no longer ac¬ 
tive and the implications of this further in Section l431 


4 DISCUSSION 

Over the last 40 years, models of spectral ageing have become 
a commonly used tool in determining the age, hence dynamics, 
of powerful FR-II galaxies. However, as was discussed in Sec¬ 
tion [T^ the validity of a number of assumptions used in the ap¬ 
plication of these models has recently been called into question. 
When more tightly constrained by modern, broad bandwidth ob¬ 


servations, |HMWood|er^ 1 I2OI3I) have shown that, for at least two 
sources, model parameters such as the injection index, and general 
assumptions such as negligible cross-lobe age variations, are less 
reliable than previously thought. These findings, combined with an 
uncertainty as to which spectral ageing model best describes emis¬ 
sion from the lobes of F R-II galaxies and the well known sp ectral vs 
dynamical age disparity l lEilekLri9^ ; lHarwood et aul2013h . means 
that if spectral ageing is to remain a valid tool for the investigation 
of such sources, it is vital that the cause of these problems are well 
understood. In this paper, we have presented results whic h further 
expand on the sample presented bv iHarwood et alj ( 1201 3h and, for 
the first time, tightly constrain curvature at frequencies of a few 
GHz using truly broad bandwidth observations. In the following 
sections we relate these new results to previous findings and dis¬ 
cuss how this changes our overall understanding of radio galaxy 
dynamics. 


4.1 Injection index 

One of the key findings of lHarwood et alj 1 I2OI3I) was the deviation 
of the injection index from previously assumed values of around 
0.5 to 0.6. It is again clear from the results presented in Section lUTI 
that this discrepancy remains, with values ranging between 0.78 
(KP) and 0.84 (JP) for 3C438 and between 1.12 (KP) and 1.21 (JP) 
for 3C28. It was noted in Section lrTI that for 3C438, at least some 
of this difference is likely to result from the inclusion of emission 
coincident with the strongly interacting jets. There is no a priori 
reason to believe this emission will be well fitted by models of 
spectral ageing and so could be forcing the derived injection in¬ 
dex to steeper values; however, from Table |5] we see that even with 
the jet removed the injection index remains > 0.7 for all models 
and so cannot alone bring these values back into agreement. 

The very steep injection index of up to 1.21 for the JP model 
of 3C28 presents an even bigger challenge to explain. The most 
straightforward physical interpretation of such a steep injection in¬ 
dex is that the jet of 3C28 terminates, or terminated, in a very weak 
shock. Assuming a simple, non-relativistic hydrodynamical shock 
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Figure 6. Two-point spectral index map created in BRATS of the southern 
lobe of 3C438 betwe en 1.4 and 8 GHz. The data used is that presented by 
iTreichel et alj j2Q0lh . reimaged at 0.42 arcsec resolution. 


in which the magnetic field of the jet fluid is wealfl the injection 
index can be related to the Mach number by telandford & Eich^ 
Il987t) 


M = 


/ 


‘^OLinj + 3 

^OLinj 1 


( 6 ) 


which for 3C28 gives a Mach number of between 1.95 (JP) and 
2.05 (KP). While this is much weaker than the shock expected 
from the historically typical injection index value of ainj = 0.6 
(M = 4.6), lower Mach numbers are n ot without precedent. For 
example, a study bv lO’Dea et al.l ( l2009l) who used high resolution 
measurements of hotspots to determine the injection index, found 
values ranging between ainj = 0.52 and 1.11 for there sample, 
corresponding to a Mach number of between 2 and 10. 

An alternative explanation is the reliability of the minimiza¬ 
tion method use d to determine the injection index. Tests on sim¬ 
ulated data (e.g. IStroe et all l2014h have shown that the min¬ 
imization used is robust and provides a significant improvement 
over previous methods such as the use of the spectral index mea¬ 
sured in particle acceleration regions; however, while this method 
is able to use all observed emission in determination of the injec¬ 
tion index, rather than just that a small zero age region, its accu¬ 
racy becomes increasingly model dependent as the fraction of aged 
plasma increases. The lack of low age emission could therefore be 
influencing the result if the model is not sufficiently constrained. 


particularly at lower frequencies. While this may provide some ex¬ 
planation of the very steep injection index observed in 3C28, from 
Figures |7] and and Table we see that the models do provide a 
good overall description of the source and are therefore unlikely to 
account for the full discrepancy. 

These findings therefore appear to reinforce the idea that 
steeper than expected injection indices are widespread. Initial find¬ 
ings from spectral studies of similar FR-II sources using low- 
frequency (50 - 160 MHz) LOFAR observations (Harwood et. al., 
in prep) also show that these steeper injection indices are observed 
at low energies where negligible losses mean that the spectral in¬ 
dex of a source should be close to the injection index. Given the 
lack of strong, distinct hotspots in our observations (and there¬ 
fore presumably weak jets), the steep injection indices we mea- 
sure are physically plausib le in models such as those presented by 
iKonar & Hardcastla l2013ll which suggest that jet power is the pri¬ 
mary factor in determining a source’s injection index. However, 
although we can be confident the injection index values observed 
are intrinsic to the source, rather than a systematic offset, the phys¬ 
ical interpretation of why such a disparity exists requires further 
investigation to determine conclusively. 


4.2 Thejetsof3C438 

As emission coincident with the location of the jets is observable 
in 3C438, the bias they may cause to the overall study of the lobe 
emission must be considered. Comparing the statistics of Table 
for 3C438, we find that for the best fitting KP model when the jet 
is included, approximately 42 per cent of the fits are rejected at > 
99 percent confidence, and without the jet only around 37 per cent 
of model fits are rejected at this level. Conversely, only 26 per cent 
of the regions show a good fit at the < 68 percent confidence level, 
compared to 28 per cent when the jet is removed. It therefore ap¬ 
pears that the inclusion of the jet may give rise to at least some 
poorly fitted regions of the source, but at a level which cannot ex¬ 
plain the rejection of the models as a whole. We discuss this model 
rejection further rejection in Section lT^ 

The observed jet emission also provides the opportunity to 
consider the dynamics of the jets themselves. As noted in Section 
13.11 brightening occurs in both the northern and southern jets co¬ 
incident with the location where it enters the lobe. One possible 
explanation of such features is therefore an interaction between the 
lobe material and the jet itself. The exact cause of such an interac¬ 
tion is unclear, but it could plausibly be due to either lobe material 
being reaccelerated at the jet-lobe boundary, or the onset of turbu¬ 
lence within the jet resulting in localised shock acceleration. The 
observed (spatial) curvature in the jet found in the southern lobe 
may possibly be evidence of this interaction; however, due to the 
presence of two low spectral age hotspots-like structures (one ac¬ 
tive and one relic, Figuresl4land [jPh . we suggest that, assuming an 
FR-II type jet, the cause is more likely due to transverse movement 
of the jet between these two regions. 

An alternative possibility is that the jets are more like those 
found in FR-Is with the lobes being formed upstream of the parti¬ 
cle acceleration regions, similar to the intermediate FR-I/II source 
Hydra A. There is some evidence to support this in the southern 
lobe where we observe a flat spectral index (Figure[3( and low-age 


® In reality, the hotspot m ay be a relativistic hydrodYnamic shock of a 
highly magnetised jet fluid ikirk et ^ 120001 : iKonar & HardcastlS . l2013h . 
but these assumptions should provide a sufficiently accurate approximation 
for our purposes. 


Note that these hotspot-like structures are also observed in Figure[6] but 
are less obvious due to the colour scale being set such as to emphasise the 
artefacts discussed in Section [T2l 
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Figure 7. Spectral ageing maps (left), corresponding maps (middle) and error maps (right) of 3C28 with 7.9 GHz flux contours. Three model fits are 
shown; JP model (top), KP model (middle) and Tribble model (bottom) using the best fitting injection index given in Table[5] 
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emission with a slight gradient away from the core (Figure|4j which 
could plausibly be an FR-I type acceleration region responsible for 
forming the lobes; however, emission from the northern jet does not 
display this form of structure and is more ’knot-like’ with is a dis¬ 
tinct separation between these regions and the low age emission at 
the end of the lobe where particle acceleration must also be occur¬ 
ring. It is therefore harder to reconcile the northern lobe with a FR-I 
type jet. It is possible that the jets are asy mmetric, as is observed in 
hybri d morphology sources (HYMORS) dGopal-Krishna & Wiital 
bOOOl) . but this would require intrinsically dissimilar jets, environ¬ 
ments and/or acceleration mechanisms to conspire to form lobes 
that are similar in terms of both morphology and the distribution of 
ages and is therefore unlikely. 

While the majority of the observed jet emission appears (at 
least in projection) to be embedded within the lobe material, the 
regions closest to the core extend well beyond the visible extent of 
the lobes and so cannot be interacting with the lobe plasma. The jets 
must therefore be interacting with their external environment, with 
the observed emission being caused either by instabilities within 
the jet and the onset of turbulence (similar to that with the lobe ma¬ 
terial discussed above), or through a more direct interaction such 
as the entrainment of thermal plasma. The morphology and spec¬ 
trum of these jets alone do not allow us to determine which of these 
two possibilities is the underlying cause of the emission but, when 
combined with the ageing analysis and source dynamics discussed 
below, does allow us to provide potentially interesting insights into 
the underlying physics off 3C438. We discuss this and its poten¬ 
tial impact on outstanding questions in the context of the sources 
dynamics further in Section l43] 


4.3 The age and dynamics of 3C438 


One of the primary concerns raised in the discussion of spectral 
ageing is the well know n discrepancy between a source’s dynami - 
cal and spectral age (e.g. lEilekLll996l;lBlundell & Rawlingslboodl) . 
Thanks to X-ray investigations by Kraft et alT i 2007 ). the environ¬ 
ment in which 3C438 reside s is well known a nd so we are able to 
constrain the dynamical age. lKraft et al.l (l2007h find that the cluster 
in which 3C438 is located has a core radius of approximately 30 
arc seconds, so that the source resides entirely within a region of 
approximately constant density. The dynami cs of the radio source 
are thus particularly simple (e.g. lFalla.ll99ll) 

1/5 




,3/5 


(7) 


where L is the length of the jet, Q is the jet power, p is the envi¬ 
ronment density, t is the time the source has been active and a is 
a dimensionless constant. We can estimate a by assuming that the 
expansion of the source is ballistic, and thus r elativistic, on scales 
smaller than the characteristic scale defined bv lFalld ( ll99 ih 


Ic = 


( 8 ) 


(here we assume a relativistic jet where the effective mass flux is 
Qjc). Then at some time to we have dL/At = c and to = Ic/c, 
which allows us to find a = 5/3. In practice the dynamics at the 
transition between the two phases will be more complicated than 
this, but the simple analysis should give an approxima tely correct 
value of a. From the observations of iKraft et al.l ( l2007h the density 
can be estimated a s p = 5.7 x 10~^^k g m“®. To estimate Q we use 
the relationship oflWillott et alj l ll999l) . normalized as described by 
iHardcastle et alj ( 12007 ). which gives Q = 1.3 x 10®® W, and we 


take L = 50 kpc. Note that Ic « 30 pc, so L ^ Ic, and the analysis 
of Falle applies. Finally we have 

/ T \ / \ 1/3 

Although the va lu e o f Q is very uncertain (see 
IHardcastle & Krausel 1201 3l) the dependence on Q is weak, 
so this is reasonably robust, and t is significantly greater than 
the observed spectral age. We have assumed the source is in 
the plane of the sky: any projection effects will make L and 
therefore t larger. A large disparity between the spectral and 
dy namical ages therefore appears to exist, similar to that observed 
bv iHarwood et al.l ( 1201 3h even when tighter constraints are placed 
on the spectral curvature at GHz frequencies by the truly broad 
bandwidth obs ervations. This age di f ferenc e is in contrast to the 
suggestion by iBlundell & Rawlings! (l2000l) who propose that it 
is possible to reconcile the dynamical age of sources with ages 
< 10 Myrs. However, for this to apply it is a requirement that the 
evolution of the magnetic field over the lifetime of the source is 
well approximated by the spectral ageing model, which may not be 
the case for 3C438. Although it is not possible to trace the history 
of the magnetic field the uncertainty introduced to a low age source 
is likely to be small compared to those in older radio lobes. While 
better approximating the magnetic field may therefore play a role 
in resolving the age disparity it is unlikely to be the only significant 
factor in resolving this problem. 

With the advance speeds determined and the external environ¬ 
ment well known, a Mach number can also be derived. The sound 
speed of the medium through which the lobes are advancing is 
given by Csound = \J 'ykTjm where 7 = 5/3 is the adiabatic 
index, m = O.bnip is the mass of the particles and T is the gas 
temperature, which in turn one can use to deri ve the Mach num- 
ber by the standard equation M = iKraft et al.l ll2007h 

find that the cluster environment in which 3C438 reside s has a gas 
temperature of 10.7 keV (region 2 of iKraft et al.Ll200^ . hence us¬ 
ing the advance speeds listed in Table [T] we find Mach numbers 
of between 3.7 [2.9] (KP) and 5.4 [4.8] (JP) for the northern lobe, 
and 3.1 [2.5] (KP) and 4.5 [4.0] (JP) for the southern lobe where 
the values in square brackets denote the values with the jet emis¬ 
sion excluded. This implies that the lobes are overpressured with 
respect to the external medium and expanding supersonically. At 
these high Mach numbers one would expect to observe a strong 
shock at the advancing edge of the lobes; however, in such a high 
temperature environment such a shock is unlikely to be observable. 
The dynamical ages provide a significantly lower advance speed of 
around 3000 km s“®, equivalent to Mach « 1.6. This means that in 
both cases the lobes are overpressured although, if the dynamical 
ages are correct, there is no expectation of a strong shock at the 
leading edge of the lobes. 

An additional way we can investigate the lobe pressures 
of 3C438 is through its energetics. From the SYNCH code of 
IHardcastle et alj h998l) (used in Section 12.4.21 to determine the 
magnetic field strength), we find the total energy density for each 
radio lobe is between 8.88 x 10“®® (KP, excluding jet) and 1.28 x 
10“®® (JP, including jet) J m“®, with a subsequent internal pres¬ 
sure, given by P = (7/3 where U is the total energy density, of 
2.96 and 4.28 x 10“®® Pa. The gas pressure of the external medium 
can also be derived by the standard equation P = nkT where 
n « 2.2no is the number density of particles. Using the gas tem¬ 
perature of 10.7 ke V and number den sity of protons of 3.4 x 10“® 
cm“® presented bv iKraft et al.l ( |2007|) . we derive the pressure of the 
external medium to be 1.3 x 10“®® Pa. These calculations therefore 
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Figure 8 . Plots of flux against frequency of example well fitted regions for 3C438 (left) and 3C28 (right), at low (top), moderate (middle) and high (bottom) 
spectral ages. Red lines indicate the best-fitting JP model with pai'ameters as noted on each plot. 


© 0000 RAS, MNRAS 000, 000-000 








16 J.J. Harwood et al. 


imply that the lobes are underpressured by over an order of mag¬ 
nitude. Such underpressures have also been measured in similar 
sources su ch as Cygnus A that are fou nd in rich cluster environ¬ 
ment (e.g. lHardcastle & Crostonl . l2010h . although the discrepancy 
is usually smaller making 3C438 the most underpressured FR-II 
source known to date. 

One possible explanati on for this difference is the uncertainty 
in the X-ray measurements, [^hellenberger et al. 1 120151) show that, 
particularly at high energies, there is a diffe rence between the tem¬ 
peratures measured by Chandra (used by iKraft et ^ l2007h and 
XMM, although it is not yet clear which of these instruments is 
correct. Assuming the lower XMM value and the offset between 
th e two instruments reduces t he temperature to 10 keV (figure 7 
of ISchellenberger et al.l ( l2015h ) we find that the pressure only re¬ 
duces to 1.2 X Pa, so remains significantly underpressured 

(we would require a temperature of only around 1 keV to bring the 
lobes back in to pressure balance). It is therefore highly unlikely 
that instrumental effects alone can resolve this issue. 

A more physical explanation for the disparity is that the lobe 
energy density is dominated by non-radiating particles. As was 
discussed in Section the most plausible cause of the strong 
emission coincident with the jets is due to an interaction with the 
external environment. Such an interaction leads to the possibility 
that thermal matter is being entrained, si milar to what is thought 
to occur in FR-I t ype radio galaxies (e.g. IPerucho & Martil200^ : 
iLaing et^l2008l) . As the jets of FR-II galaxies must remain highly 
relativistic out to l arge distance s, entrainment is thought to be low 
in most cases (e.g. iLeahvl 1 199 ih : however, 3C438 does not display 
the strong, compact hotspots that are typical of FR-II galaxies. If 
for 3C438 the cause of these weak hotspots is due to a relatively 
slow moving jet, then conditions would be favourable for entrain¬ 
ment to occur. In such a scenario, the assumed minimum energy 
condition would no longer apply, significantly altering the spectral 
ages and pressures derived for the source. Such pressure dispar¬ 
ities are also observed in other sources where material is plausi- 
bly being entrained (e . g. the intermediate FR-I/II source Hydra A, 
iHardcastle & CrostonL 201(]ti and so may be relatively common in 
lobes sources with weak jets. A slow moving jet is also able to pro¬ 
vide a natural solution with respect to the observed initial electron 
energy distributio n as, assuming r, et < 10, a steep injection index 
is to be expected jKirk et alil200(]ll^^L 

Another possible explanation comes from investigations at X- 
ray energies which have shown that the magnetic field strength of 
many F R-II sources is not i n equipartition as has been so far as¬ 
sumed dCroston et all l2005h . Values range from between 0.3 and 


ing our determination of the internal lobe pressure to be underesti¬ 
mated. Using the IP model with the jet included (where the internal 
pressure is highest) a nd assuming the minimum value of 0.3 of 
ICroston et al. l l2005h . we find that the lobes remain underpressured 
by an order of magnitude, with the pressure increasing to only 1.98 
xl0““ Pa. 

We also know that a maximum spectral age occurs when B = 
Bcmb/'/S llHarwood et '^l2013h where Br-Mn = 0.31 8(14-^:)^ 
nT is the equivalent magnetic field strength of the CMB l iHughesl . 
Il99lh which, at the redshift of 3C438, occurs when B = 0.31 Beq. 
This is close to the value used to derive the maximum pressure and 


Note that this does not explicitly exclude the poss ibility of a fast jet for 
which a steep injection index can also be produced jKonar & HardcastlS . 
| 2013 ^ 


increases the spectral age of 3C438 to «19 Myrs. Therefore, while 
departure from equipartition provides a possible explanation for the 
disparity between the spectral and dynamical ages, it is unlikely 
alone to provide the solution to the underpressured lobes of 3C438. 
Determining the cause of such disparities will be vital for future 
studies if we are to understand the dynamics and energetics of these 
sources. 


4.4 The classification of 3C28 

While 3C438 presents a morphology typical of many archetypal 
FR-IIs, 3C28 has historically been much more difficult to classify. 
While currently classified in the 3CRR catalogue as an FR-II, its 
non-standard morphology and lack of strong hotspots mean that in- 


vestigations have both historically (e.g.lMacklin 

.ll983l;lLaing et all 

1 19831; Hardcastle et al.L 19991; Chiaberge et all 

1999h and recently 


(e.g. IPonato et all |2004 : iBalmaverde et al.L 120061) deemed 3C28 


to be an FR-I; a classification which also persists in major astro¬ 
nomical databases (e.g. the NASA Extragalactic DatabasJ^. As 
within this paper we present both the most sensitive, high resolu¬ 
tion images at GHz frequencies of 3C28 to date, along with for the 
first time detailed information about the source’s spectrum on small 
spatial scales, we briefly consider which of these classifications is 
correct. 

While at low resolutions the unusual morphology of 3C28 un¬ 
derstandably causes problems in terms of classification, from Fig- 
ure[T]we see that when well resolved the source’s overall morphol¬ 
ogy closely resembles that of an FR-II. The double lobe structure 
displays significant edge brightening with only a thin jet or bridge 
observable near the centre of the source and although the northern 
lobe lacks any compact structure, the southern lobe contains what 
closely resembles a hotspot, albeit a weak one. The spatial distri¬ 
bution of spectral age shown in Figure |7] also supports the FR-II 
case, with the gradient of ages increasing towards the core being 
indicative that the particles were accelerated at the tip of the lobe, 
rather than close to the core region as is observed in FR-I sources. 
We are therefore confident that the FR-II classification which has 
been assumed throughout this paper is robust. 


4.5 Where in the radio lifecycle is 3C28? 

While as discussed in the previous section the distribution of spec¬ 
tral ages are as expected for an FR-II galaxy, they also raise an 
interesting question; where exactly in the radio lifecycle is 3C28? 
Intuitively, one would expect the distribution to be arranged such 
that zero or low age emission is observed in or close regions of par¬ 
ticle acceleration which then ages with distance from the hotspot, 
as is seen in 3C438 (Figure S and in those sources studied by 
iHarwood et al.l i2013h . However, for 3C28 this is not the case and 
we observe the minimum age to be on the order of around 6 to 9 
Myrs (Figure|7j. One possible explanation for this is the combina¬ 
tion of a relatively weak shock and the unavoidable consequence of 
measuring emission from a range of plasma ages along the line of 
sight. This superposition of spectra is usually assumed to be negli¬ 
gible but if emission along the line of sight to the hotspot is particu¬ 
larly significant it may dominate and cause no low age regions to be 
observed. While this provides a plausible explanation for the north¬ 
ern lobe where all of the observed emission is highly diffuse, the 
presence of the bright compact structure in the southern lobe makes 


https://ned.ipac.caltech.edu/ 
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this an unlikely scenario. Emission from the southern hotspot dom¬ 
inates the total flux of the region and therefore, assuming current 
particle acceleration observed in superposition with older regions 
of plasma, should still display a reduced age relative to the sur¬ 
rounding plasma, even if that age is non-zero. The smooth distri¬ 
bution of spectral ages across both diffuse and compact structure 
therefore suggests an alternative scenario. 

While the ‘jet-like’ emission close to the centre of the source 
has previously led to the belief that 3C28 is currently active, the 
lack of a radio (Figure [TJ or X-ray core jHardcastle et all 1 20061) 
supports the idea that the central AGN is turned off, which provides 
a natural explanation for the observed smooth spatial distribution of 
spectral ages. Once the central engine becomes inactive and freshly 
accelerated particles are no longer supplied via the hotspots, the 
overall brightness of the lobes will begin to fade due to both ra¬ 
diative and adiabatic loses. As adiabatic losses are independent of 
frequency only the brightness of the observed emission is affected 
and not the curvature within its spectrum. If therefore the adia¬ 
batic losses remain low in the region of the (now inactive) hotspot, 
one will observe relatively bright, compact structure that is highly 
curved such as is found in the southern lobe of 3C28. We there¬ 
fore suggest that 3C28 is currently inactive with the AGN having 
shut down between 6 and 9 Myrs ago, and the ‘jet-like’ emission 
close to the centre of the source being a result of a bridge similar 
to those observed in the relaxe d double 3C442A dComins & OwenL 
ll99lHHardcastle et ahlflg^ . or jet interaction during a previous 
episode of AGN activity. 


4.6 Model comparison 

Determining which, if any, model of spectral ageing best describes 
the observed emission from FR-IIs is a key step in producing reli¬ 
able ages for these sources. The model fits of 3C28 provide good 
agreement with the observed spectrum with none of the models 
being rejected at the 68 per cent confidence level (Table [^; how¬ 
ever, the goodness-of-fit for 3C438 fares much worse. Given the 
fairly typical morphology and spectral age distribution of 3C438, 
one would expect the three models of spectral ageing tested to be 
well fitted to the observations but are all instead rejected at the 99 
per cent confidence level. As was discussed in Section |42] emis¬ 
sion coincident with the jet may cause some bias towards higher 
values but is unable to fully account for the model rejections. 

Evidence of another cause of this poor goodness-of-fit lies in 
the northern hotspot region, where strong emission is observed but 
the models are poorly fitted. One potential cause of this is a dif¬ 
fering injection index in the hotsp ot and the lobes, as has been ob- 
served in other FR-II sources (e.g. lHardcastle & LoonevLl200lh . If 
particle acceleration occurs over an extended region rather than at 
a single point, the lobe emission may be best described by the elec¬ 
tron distribution as it leaves the hotspot, rather than in the hotspot 
itself where turbulence and the mixing of electron populations is 
likely to be high. The hotspots in the southern lobes do not display 
such a reduction in the goodness-of-fit, but given this emission is 
much weaker than its northern counterpart, this is perhaps not sur¬ 
prising. High resolution studies designed to better understand the 
workings of FR-II hotspots may therefore be required if one wishes 
to eventually determine a framework for highly reliable, detailed 
models of spectral ageing. 

Although jet and hotspot emission are likely to cause some 
bias to high values, insight into the underlying reason for 
the overall model rejection can also be found by comparing the 
goodness-of-fit between the two lobes of 3C438. The extended 


emission in the northern lobe is generally well described by the 
spectral ageing models compared to its southern counterpart which 
overall has much higher values. The worst fitting regions of ex¬ 
tended emission in the southern lobe appear to be located between 
the two hotspots and close to the path of the jet. As the jet and 
hotspot sweep across the lobe between the two areas, freshly accel¬ 
erated particles will become mixed with older regions of plasma. 
This mixed electron population will not be well described by the 
single injection models used within this paper and we therefore 
suggest this to be the primary cause of poor model fits, rather than 
any inherent problem with the models themselves. Detailed numer¬ 
ical modelling of spectral ages in radio galaxy lobes may be able to 
reproduce such effects in the future. 

An additional point which must be addressed is; which of the 
three models tested provides the best description of these sources 
and is such a model physically realistic? From Table|^we see that 
for both 3C438 and 3C28, the KP model provides the best descrip¬ 
tion of the observed spectra, but an environment in which the pitch 
angle remains fixed (the characteristic feature of the KP model) 
requires a set of finely tuned parameters and therefore provides 
a less likely physical interpretation. Whilst historically sources in 
which t he KP mode l provi des the best description are not unheard 
of (e.g. lCarilli et alill99lh . they were in general thought to be the 
exception. However, we find that when considered on small spatial 
sc ales, all of the FR-II s ources tested both within this paper and 
bv iHarwood et alj ) l20I3h show that the KP models provides a sys¬ 
tematically better fit than the more physically plausible IP model, 
suggesting this may be common in the majority of sources. Deter¬ 
mining whether these results are due to intrinsically fixed pitch an¬ 
gles, or if improved time-averaged pitch angle models are required, 
is key to the ability of spectral ageing in providing a more accu¬ 
rate description of the spectra of FR-II sources in the future. The 
Tribble model is able to provide (at least to some extent) an an¬ 
swer to this problem, giving the physical realism of the IP model, 
an improved goodness-of-fit and a more advanced description of 
the magnetic fields within the lobes. However, it is clear from these 
results that if we are to begin eliminating spectral ageing models 
at a statistically significant level, one must be both well sampled in 
frequency space, so as to tightly constrain model parameters, and 
fitting performed over very wide frequency range where the spec¬ 
tral variations between the 3 models are significant. 


5 CONCLUSIONS 


In this paper we have presented a spectral ageing study of two 
cluster-centre FR-II sources, 3C438 and 3C28. Using broad band¬ 
width VLA observations to tightly constrain spectral curvature at 
GHz frequencies, we have investigated their dynamics and ener¬ 
getics and determined which model of spectral ageing provides the 
best description of these sources. We have looked at whether sig¬ 
nificant cross-lobe variations are likely to be common in powerful 
radio source s and if the high inje ction indices found in our pre¬ 
vious work jHarwood et al. . |20 1 3l) are still observed when tightly 
constrained by full bandwidth observations. We have investigated 
whether emission from the jet of 3C438 affects spectral ageing 
analysis and test whether the well known dynamical versus spec¬ 
tral age disparity still exist in cluster-core FR-IIs. The key points 
made within this paper are as follows: 


(i) We find that the injection index of both sources are higher 
than previously assumed, consistent with those found in our 
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previous investigations. 

(ii) Jet emission causes a bias towards steeper injection index 
values but, even for the strong Jet emission observed in 3C438, 
does not account for large difference between the commonly 
assumed values and those found within this paper. 

(iii) Significant cross-lobe age variations are again observed in 
the lohes of 3C438. We suggest that these variations are a common 
feature of FR-II sources and so should be considered carefully in 
all future studies. 

(iv) We suggest that Tribble model currently provides the most 
convincing description when both goodness-of-fit and physical 
plausibility are considered for both 3C438 and 3C28. 

(v) We derive a Mach number for the lohes of 3C438 of 
between 3.1 and 5.4 (1.6 for the dynamical age), implying that 
the radio lobes are overpressured with respect to the external 
medium. However, the source energetics suggest that the radiating 
particles and magnetic field at equipartition cannot account for the 
necessary pressure to support the lobes, similar to what is observed 
in other rich cluster sources (e.g. Cygnus A). A discrepancy of 
approximately an order of magnitude makes 3C438 the most 
underpressured FR-II source known to date. 

(vi) We find that even when tightly constrained at GHz fre¬ 
quencies, an order of magnitude disparity between the spectral and 
dynamical ages remains. We show that for 3C438, both departure 
from equipartition and the presence of non-radiating particles 
arising from the entrainment of thermal material are both able to 
provide a possible solution to this problem. 

(vii) We confirm that 3C28 is an FR-II, rather than an FR-I as 
reported in some previous investigations, based on its morphology 
and spectrum at high resolution. 

(viii) We suggest that 3C28 is a relic source with the central 
AGN turning off around 6 to 9 Myrs ago. 

We therefore confirm that many of the assumptions which 
are made in relation to spectral ageing models should be carefully 
reevaluated. The underlying physical reason for these disparities, 
such as the high injection indices, remains an open question which 
will be addressed in the near future by low-frequency studies of ra¬ 
dio lobes and of hotspots in FR-IIs (e.g. Harwood et al., in prep). 
We conclude that if one is to determine reliable intrinsic ages, and 
hence the total power of these sources, then it is vital that these out¬ 
standing questions are resolved and, possibly, new models which 
better account for our expanding knowledge of these sources de¬ 
rived. 


6 ACKNOWLEDGEMENTS 

We wish to thank the anonymous referee whose constructive sug¬ 
gestions have helped expand our interpretation of the results. JJH 
wishes to thank the Netherlands Institute for Radio Astronomy 
(ASTRON) for a postdoctoral fellowship and the University of 
Hertfordshire and the Science and Technology Facilities Coun¬ 
cil (STFC) for funding via the Studentship Enhancement Pro¬ 
gramme (STEP) award. JHC is grateful for support from the Sci¬ 
ence and Technology Eacilities Council under grant ST/JOO1600/1. 


This work has made use of the University of Hertfordshire Sci¬ 
ence and Technology Research Institute high-performance comput¬ 
ing facility. This research has made use of the NASA/IPAC Extra- 
galactic Database (NED) which is operated by the Jet Propulsion 
Laboratory, California Institute of Technology, under contract with 
the National Aeronautics and Space Administration. We wish to 
thank staff of the NRAO Jansky Very Large Array of which this 
work makes heavy use. The National Radio Astronomy Observa¬ 
tory is a facility of the National Science Foundation operated under 
cooperative agreement by Associated Universities, Inc. 


References 

Alexander P, 1987, MNRAS, 225, 27 
Alexander P, Leahy J. P, 1987, MNRAS, 225, 1 
Avni Y., 1976, ApJ, 210, 642 

Balmaverde B., Capetti A., Grandi P, 2006, A&A, 451, 35 
Beers T. C., Geller M. J., Huchra J. P, 1983, ApJ, 264, 356 
Begelman M. C., Cioffi D. F, 1989, ApJ, 345, L21 
Birkinshaw M., Laing R., Peacock J. A., 1981, MNRAS, 197, 253 
Blandford R. D., Eicher D., 1987, PhR, 154, 1 
Blandford R. D., Ostriker J. P, 1978, ApJ, 221, 29 
Blundell K. M., Rawlings S., 2000, ApJ, 119, 1111 
Blundell K. M., Rawlings S., 2001, in Laing R., Blundell K., eds, 
ASP Conference Series Vol. 250, ‘Particles and Field in Radio 
Galaxies’. 

Bower R. G., Benson A. J., Malbon J. C., Helly J. C., Erenk C. S., 
Baugh C. M., Cole S., Lacey C. G., 2006, MNRAS, 370, 645 
Burch S. F, 1977, MNRAS, 180, 623 
Burch S. F, 1979, MNRAS, 186, 519 
Carilli C., Perley R., Dreher J., Leahy J., 1991, ApJ, 383, 554 
Chiaberge M., Capetti A., Celotti A., 1999, A&A, 349, 77 
Comins N. R, Owen F. N., 1991, ApJ, 382, 108 
Cornwell T. J., 2008, JSTSP, 2, 793 

Croston J. H., Hardcastle M. J., Harris D. R., Belsole E., Birkin¬ 
shaw M., Worrall D. M., 2005, ApJ, 626, 733 
Croton D. J. et al., 2006, MNRAS, 365, 11 
de Gasperin F. et al., 2012, A&A, 547, A56 
Donato D., Sambruna R. M., Gliozzi M., 2004, ApJ, 617, 915 
Eilek J. A., 1996, in Hardee P. E., Bridle A. H., Zensus J. A., 
eds, ASP Conference Series Vol. 100, ‘Energy transport in radio 
galaxies and quasars’. 

Eilek J. A., Arendt P. N., 1996, ApJ, 457, 150 
Eilek J. A., Melrose D. B., Walker M. A., 1997, ApJ, 483, 282 
Falle S. A. E. G., 1991, MNRAS, 250, 581 
Fanaroff B. L., Riley J. M., 1974, MNRAS, 167, 31P 
Gopal-Krishna Wiita P. J., 2000, A&A, 363, 507 
Gregorini L., Bondi M., 1989, A&A, 225, 333 
Hardcastle M. J., 2013, MNRAS, 433, 3364 
Hardcastle M. J., Alexander P, Pooley G. G., Riley J. M., 1997, 
MNRAS, 288, 859 

Hardcastle M. J., Alexander P, Pooley G. G., Riley J. M., 1999, 
MNRAS, 304, 135 

Hardcastle M. J., Birkinshaw M., Worrall D. M., 1998, MNRAS, 
294, 615 

Hardcastle M. J., Croston J. H., 2010, MNRAS, 404, 2018 
Hardcastle M. J., Evans D. A., Croston J. H., 2006, MNRAS, 370, 
1893 

Hardcastle M. J., Evans D. A., Croston J. H., 2007, MNRAS, 376, 
1849 


© 0000 RAS, MNRAS 000, 000-000 


Spectral ageing in cluster-centre FR-lIs 19 


Hardcastle M. J., Kraft R. R, Worrall D. M., Croston J. H., Evans 

D. A., Birkinshaw M., Murray S. S., 1997, MNRAS, 662, 166 
Hardcastle M. J., Krause M. G. H., 2013, MNRAS, 430, 174 
Hardcastle M. J., Looney L. W., 2001, MNRAS, 320, 355 
Hardcastle M. J., Massaro H., Harris D. E., 2010, MNRAS, 401, 

2697 

Hardcastle M. J., Sakelliou 1., 2004, MNRAS, 349, 560 
Harwood J. J., Hardcastle M. J., Croston J. H., Goodger J. L., 
2013, MNRAS, 435, 3353 

Heesen V., Croston J. H., Harwood J. J., Hardcastle M. J., Ananda 
H., 2014, MNRAS, 439, 1364 

Hughes P. A., 1991, Beams and Jets in Astrophysics. Cambridge 
University Press 

Jaffe W., Perola G., 1973, A&A, 26, 423 

Jamrozy M., Machalski J., Mack K.-H., Klein U., 2005, A&A, 
433, 467 

Kaiser C. R., Alexander R, 1997, MNRAS, 286, 215 
Kardashev N. S., 1962, AJ, 6, 317 

Kharb R, O’Dea C. R, Baum S. A., Daly R. A., Mory M. R, Don¬ 
ahue M., Guerra E. J., 2008, ApJS, 174, 74 
Kirk J. G., Guthmann A. W., Gallant Y. A., Achterberg A., 2000, 
ApJ, 542, 235 

Komissarov S. S., Gubanov A. G., 1994, A&A, 285, 27 
Konar C., Hardcastle M. J., 2013, MNRAS, 436, 1595 
Konar C., Saikai D. J., Jamrozy M., Machalski J., 2006, MNRAS, 
372, 693 

Kraft R. R, Eorman W. R., Hardcastle M. J., Jones C., Nulsen P. 

E. J., 2007, ApJ, 664, 83 

Krause M., Alexander R, Riley J., Hopton D., 2012, MNRAS, 
427, 3196 

Laing R. A., Bridle A. H., Cotton W. D., Worrall D. M., Birkin¬ 
shaw D. M., 2008, in Rector T. A., De Young D. S., eds, ASP 
Conference Series Vol. 386, ‘Extmgalactic jets: theory and ob- 
serx’ations from radio to gamma ray’. 

Laing R. A., Riley J. M., Longair M. S., 1983, MNRAS, 204, 151 
Leahy J. P., 1991, Beams and jets in astrophysics: Interpretation 
of large scale extragalactic jets. Cambridge University Press 
Leahy J. R, Black A. R. S., Dennett-Thorpe J., Hardcastle M. J., 
Komissarov S., Rerley R. A., Riley J. M., Scheuer R. A. G., 1997, 
MNRAS, 291, 20 

Machalski J., Jamrozy M., Saikia D. J., 2009, MNRAS, 395, 812 
Macklin J. T., 1983, MNRAS, 203, 147 
Marscher A. R, Jorstad S. G., 2011, ApJ, 729, 26 
Meisenheimer K., Rbser H. J., Hiltner R. R., Yates M. G., Longair 
M. S., Chini R., Rerley R. A., 1989, A&A, 219, 63 
Mullin L. M., Hardcastle M. J., Riley J. M., 2006, MNRAS, 372, 
113 

O’Dea C. R, Daly R. A., Kharb R, Ereeman K. A., Baum S. A., 
2009, A&A, 494, 471 

O’Dea C. R, Owen E. N., 1986, ApJ, 301, 841 
O’Dea C. R, Sarazin C. L., Owen E. N., 1987, ApJ, 316, 113 
Orru E., Murgia M., Eeretti L., Govoni R, Giovannini G., Lane 
W., Kassim N., Raladino R., 2010, A&A, 515, A50 
Owen E. N., Ledlow M. J., 1994, in Bicknell G. V., Dopita M. A., 
Quinn P. J., eds, ASP Conference Series Vol. 54, ‘The first 
Stromlo symposium: The physics of active galaxies’. 

Owen E. N., Rudnick L., 1976, ApJ, 205, LI 
Pacholczyk A. G., 1970, Radio astrophysics. Nonthermal pro¬ 
cesses in galactic and extragalactic sources. San Francisco, Free¬ 
man 

Perley R. A., Rbser H. J., Meisenheimer K., 1997, A&A, 328, 12 
Perucho M., Marti J. M., 2007, MNRAS, 382, 526 



00S6""05= 04'' 03" 02" 

J2000 Right Ascension 


Figure Al. Combined B- and C-configuration radio map of 
J00560226+2627287 between 4 and 8 GHz. Imaged using multiscale 
CLEAN and CASA nterms = 2 to a central frequency of 6 GHz. The 
off-source RMS of the combined map is 4/rJy beam~^. The restoring beam 
is 0.99 arcsec. 


Rau U., Cornwell T. J., 2011, A&A, 532, 71 
Schellenberger G., Reiprich T. H., Lovisari L., Nevalainen J., 
David L., 2015, A&A, 575, 30 
Scheuer P. A. G., 1974, MNRAS, 166, 513 
Skrutskie M. F. et al., 2006, AJ, 131, 1163 
Spergel D. N. et al., 2003, ApJS, 148, 175 
Shoe A., Harwood J. J., Hardcastle M. J., Rbttgering H. J. A., 
2014, MNRAS, 445, 1213 

Treichel K., Rudnick L., Hardcastle M. J., Leahy J. R, 2001, ApJ, 
561,691 

Tribble R, 1993, MNRAS, 261, 57 

Willott C. J., Rawlings S., Blundell K. M., Lacy M., 1999, MN¬ 
RAS, 309, 1017 

Winter A. J. B. et al., 1980, MNRAS, 192, 931 
Zabludoff A. L, Huchra J. R, Geller M. J., 1990, ApJS, 74, 1 


APPENDIX A: J00560226+2627287: A TAILED SOURCE IN 
T HE EIELD OF 3C28 

Within the field of view of our observations another large, bright, 
extended radio source located ~4 arcminutes (~750 kpc projected) 
north east of 3C28 was also serendipitously observed (Figure |aT]i. 
Unfortunately due to being located at the edge of the JVLA’s pri¬ 
mary beam and with no well matched L-band data available, we are 
unable to perform a full spectral analysis; however, we are still able 
to provide some insight into the source’s properties. As this source 
is located well away from the phase centre, in order to determine 
reliable flux values a correction for primary beam attenuation was 
required. Each image was therefore corrected using CASA’s IMPB- 
COR task, ^^s for the JVLjA BprimaTy 45 /a (where Bp.f{. 7 nary 
is in arc minutes and v is in GHz), at higher frequencies the source 
is located very close to the edge of the primary beam with some 
regions falling outside of the corrected image. We therefore restrict 
our analysis to a frequency range of between 4 and 6 GHz to ensure 
all of the source is fully encompassed. Integrated flux values were 
obtained using CASA from which we derive a spectral index be- 
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Table Al. Integrated flux for 1005603+262717 


Frequency 

(GHz) 

RMS 

(Miy) 

Integrated Flux 
(mJy) 

± 

(mly) 

4.33 

33.6 

95.4 

1.0 

4.45 

36.3 

93.2 

0.9 

4.58 

37.3 

91.9 

0.9 

4.71 

32.8 

89.8 

0.9 

4.84 

35.5 

87.7 

0.9 

4.97 

34.5 

84.9 

0.9 

5.07 

38.9 

83.9 

0.8 

5.20 

33.4 

81.3 

0.8 

5.33 

38.9 

78.8 

0.8 

5.45 

38.1 

77.5 

0.8 

5.58 

39.7 

74.7 

0.7 

5.71 

36.2 

72.8 

0.7 

5.84 

38.9 

71.5 

0.7 

5.97 

43.3 

71.0 

0.7 


Integrated fluxes for 1005603+262717 between 4 and 6 GHz. The RMS 
noise is taken from the image after primary beam correction, in a region 
away from but representative of the source. The spectral index as determined 
by a least-squares fit to the values listed above is o = 0.98. 

tween 4.3 and 6.0 GHz of a = 0.98 using a standard least-squares 
fit. The integrated flux values and their associated uncertainties are 
shown in Table lAl] 

The host galaxy of this source has previously been identified 
as 2MASX J00560226+2627287 {z = 0.191399) and is associ¬ 
ated with_Aben ll^the_same_clusto ofwhich 3C28 is a mem¬ 
ber teeers et alill983l : IZabludoff et al.Lll990l) . From Figure lATI we 
see that the source contains four major components; a compact 
core in the north west, a pair of jet-like structures which termi¬ 
nate at bright lobes, and two diffuse tails which, in the plane of 
the sky, run parallel to the lobes of 3C28. At first sight, the mor- 
pho logy resembles that of a narrow angle tailed radio galaxy (NAT; 
e.g. lO Pea & OwenL 1 19861 : lO’Dea et all Il987h . NATs are nearly 
always found in rich clusters and thought to be a result of ram pres¬ 
sure bending the jet and lobes as they move through the intraclus¬ 
ter medium (ICM), forming their characteristic narrow tails. This 
was also the conclusion reached in a previous low -resolution radio 
study of this source bv iGregorini & Bondil ( Il989h ; however, at the 
end of the northern tail we also see diffuse emission perpendicu¬ 
lar to the main structures which is hard to explain in the context 
of a NAT as, in the rest frame of the source, it would lie in the 
flow of material in the external med i um. T his diffuse feature was 
also observed bv iGregorini & Bondil ( Il989ll and so we can be con¬ 
fident that this emission is real, rather than an imaging artefact. The 
most plausible explanation is therefore t hat the source is in s tead a 
wide-angle tail radio galaxy (WAT; e.g. lOwen & Rudnickl Il976l : 
iHardcastle & SakelliouLlrOOdh seen in projection. 

Similar to NATs, WATs are also observed to reside in rich 
cluster environments close to the core region, but are thought to 
be moving more slowly through the I CM resulting in a less pro- 
nounced bending of their jets and lobes. I&egorini & Bondil l ll989h 
derive a velocity of the host galaxy with respect to the mean cluster 
velocity of 2400 km s“^, much faster than the velocity of ~100 
km s“^ expected for WATs; however, the minimum physical sepa¬ 
ration between this source and 3C28 (located at the centre of Abell 
115) is ~750 kpc (this is really a lower limit as the redshifts sug¬ 
gests the potential WAT is significantly closer than 3C28). This far 
away from the centre, the source almost certainly resides in a lower 


density environment than those WATS observed close to the cluster 
core, which would allow jets to remain only moderately distorted 
at these high velocities. 

Another characteristic feature of WATs is that they have so 
far all been observed to be hosted by massive elliptical galaxies. 
From measurements m ade by the 2MASS survey of the host galaxy 
dSkrutskie et al.L l2006h . we see that both the potential WAT and 
3C28 (which is known to reside in a massive cD galaxy) have K- 
band magnitudes of around 13.5, suggesting that the WAT’s host 
galaxy is also similarly massive. This suggests that while WATs 
have a preference to form in cluster cores due to the abundance of 
massive elliptical galaxies, it is not a strict requirement. Further in¬ 
vestigation is required to substantiate this claim but it provides at 
least some tentative evidence that WATs in fast moving, massive 
galaxies further away from the core are also possible given the cor¬ 
rect environment. 

One further piece of evidence that supports the WAT classifi¬ 
cation comes from our radio luminosity measurements. WATs are 
thought to be a subtype of FR-II or transitional FR-I/II radio galax¬ 
ies. Assuming the spectral index given above and the flux measure¬ 
ments of Table I aTI we find a 178 MHz l uminosity of 2 x 1 0^^ W 
Hz”’^, which lies above the FR-I/II break dOwen & Ledlowl [19941) . 
With the northern lobe angled towards the observer, one can also 
naturally explain the observed sharp turn in the diffuse emission at 
the end of the northern lobe, the asymmetry between the two jet¬ 
like structures, and the northern jet’s increased flux density close 
to the core (due to beaming) through simple geometry. We there¬ 
fore conclude that the source is most likely a WAT hosted by fast 
moving, massive elliptical galaxy in the cluster Abell 115. 
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